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The self-similarity hypothesis claims that in classical general relativity, spherically symmetric 
solutions may naturally evolve to a self-similar form in certain circumstances. In this context, the 
validity of the corresponding hypothesis in nonspherical geometry is very interesting as there may 
exist gravitational waves. We investigate self-similar vacuum solutions to the Einstein equation 
in the so-called whole-cylinder symmetry. We find that those solutions are reduced to part of the 
Minkowski spacetime with a regular or conically singular axis and with trivial or nontrivial topology 
^\ | if the homothetic vector is orthogonal to the cylinders of symmetry. These solutions are analogous 

to the Milne universe, but only in the direction parallel to the axis. Using these solutions, we discuss 

■ the nonuniqueness (and nonvanishing nature) of C energy and the existence of a cylindrical trapping 
' horizon in Minkowski spacetime. Then, as we generalize the analysis, we find a two-parameter family 
, of self-similar vacuum solutions, where the homothetic vector is not orthogonal to the cylinders in 

q ■ general. The family includes the Minkowski, the Kasner and the cylindrical Milne solutions. The 

y-^ ' obtained solutions describe the interior to the exploding (imploding) shell of gravitational waves or 

the collapse (explosion) of gravitational waves involving singularities from nonsingular initial data 
£SJ ■ in general. Since recent numerical simulations strongly suggest that one of these solutions may 

£N| ' describe the asymptotic behavior of gravitational waves from the collapse of a dust cylinder, this 

means that the self-similarity hypothesis is naturally generalized to cylindrical symmetry. 

cr 

, ^ I. INTRODUCTION 

In studying nonspherical vacuum gravitational fields, cylindrically symmetric systems have the advantages of being 

, essentially 1+1 dimensional, and, unlike the spherical case, possessing a dynamical degree of freedom in gravity, 

' i.e., gravitational waves. Solutions in this system are discussed in [l[ and are called Einstein- Rosen waves. These 

, have played an important part in the history of gravitational wave research, principally in elucidating the reality of 

(*C) • gravitational waves as carriers of energy Q ■ Thus, several researchers have studied vacuum and nonvacuum cylindrical 

\ systems in an attempt to clarify the nature of (especially) nonspherical gravitational collapse, in a system that involves 

I ■ the essential nonlinearity of the gravitational field and the emission of gravitational waves but requires the anal ysis 

QO ' °f partial differential equations with just one spatial dimension. See for example pj 13. Rj FtL M, \$L fill [Til. H3, [l3l.Tl3. 

o: Eosoa. 

As in many other fundamental theories, self-similarity plays an important role in gravitation. This importance is 
• i-h , encapsulated in Carr's self-similarity hypothesis [H, [lj|, which originally asserts that in the cosmological context, 

■ spherically symmetric solutions of Einstein's equations evolve to a self-similar form. This evolution can be towards 
\ either an intermediate or an endstate. Later it was found that this is actually the case in gravitational collapse. 

Examples include the collapse of a soft fluid sphere [2(| Hl[ and critical phenomena emerging in the spherical collapse 
of a variety of matter fields (22J. We should also note that spatially homogeneous cosmological models provide an 
example of tendency towards self-similar solutions in non spherically symmetric systems: see [23j for a review. 

Spherically symmetric self-similar spacetimes have been extensively studied and are now well understood. A natural 
next step therefore is to consider the role of self-similarity in cylindrical systems. These arise as a special case of 
G2 spacetimes - i.e. spacetimes admitting a two-dimensional group of isometries, usually but not always Abelian. 
There have been several studies of self-similar G2 cosmological models - again see [23[ , but there has been little work 
heretofore on the role of self-similarity in cylindrical collapse. Indeed it is generally true that in the cylindrical case, 
results are fewer and farther between: the additional degrees of freedom that exist in the cylindrical case, which 
render the systems of equations encountered considerably more difficult, have not to date allowed the development 
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of as clear a picture as we have in the spherical case (but see [HI, [24], [2f| [26| ) . We note that Ref. [26[ deals with 
self-similar cylindrical dust collapse. 

The present work seeks to increase our understanding of the role of self-similarity in cylindrical gravitational collapse, 
and the veracity of the similarity hypothesis in this context. In particular, this work is motivated by numerical results 
obtained by two of the present authors (KN and TH) and their collaborators which show that self-similar cylindrically 
symmetric vacuum solutions appear to describe the asymptotic behavior of the (numerically obtained) gravitational 
field outside a collapsing hollow dust cylinder [27l |. We present here a two-parameter family of self-similar vacuum 
solutions in cylindrical symmetry. 

This paper is organized as follows. In Sec. II, we present formulations for vacuum spacetimes and self-similar 
spacetimes in the so-called whole- cylinder symmetry independently. In Sec. Ill, we solve the Einstein field equations 
for vacuum self-similar spacetimes and obtain exact solutions. The motivation for this is to twofold. First, we wish 
to study self-similar cylindrical collapse of pure gravitational waves, as representing the simplest example of self- 
similar cylindrical collapse. The hope is that in some sense, this would provide the standard model for self-similar 
cylindrical collapse in the same way that the Schwarzschild solution provides the standard model of the (endstate 
of) spherically symmetric collapse. The second motivation is to determine the description of possible exteriors of 
collapsing self-similar, cylindrical matter. We recall also the overarching motivation for this study: to determine 
self-similar solutions that may act as (intermediate) attractors of more general cylindrically symmetric spacetimes. 
In fact we show that Einstein-Rosen waves with a homothetic vector orthogonal to the cylinders of symmetry are flat 
and reduce to part of the Minkowski spacetime with or without a conical singularity and with or without nontrivial 
topology. Moreover j_we show that the above metric form of part of the Minkowski spacetime implies nontrivial (i.e. 
nonzero) C energy Q and a trapping horizon. In Sec. IV, we naturally extend the solutions to more general class 
and find that these solutions arc also self-similar where the homothetic vector is not orthogonal to the cylinders 
of symmetry in general. We consider the analytical extension of the solutions, analyze the global structure of the 
obtained spacetimes and find that these describe interesting nonlinear dynamics of gravitational waves. We discuss the 
physical interpretation of these solutions. In Sec. V, we summarize the paper. We use the units, in which G = c = 1. 



II. CYLINDRICAL SYMMETRY AND SELF-SIMILARITY 



A. Spacetimes in whole-cylinder symmetry 

For cylindrically symmetric spacetimes, we assume that there are two commuting, spatial Killing vectors (£(i),£(2)) 
such that the orthogonal space is integrable and the Killing coordinate ip, where = d/d(p, is identified at and 
2ir. We call and £( 2 ) azimuthal and translational Killing vectors, respectively. Here we shall additionallyassume 
that each of the two Killing vectors be hypersurface orthogonal, which is called whole-cylinder symmetry or the 
polarized case (To| . The circumferential radius p, the specific length I and the arcal radius r are defined as 

P 2 :=£(!)«£(■), ^ := £(2)^2) and r := pi. (2.1) 

Note then that r is the areal radius of the orbits of the isometry group and so r > with r = at the axis. The line 
element in this class of spacetimes is given by [28| : 

ds 2 = -2e 2 ~> {u < v) - 2 ^ u ' v) dudv + e- 2 ^ u ' v) r 2 {u,v)d(p 2 + e 2 ^ u ' v) dz 2 . (2.2) 

We note that this form of the line element is unchanged under rescalings of the null coordinates u — > u(u), v — > v(v). 
The regular axis condition fiol . [28j ensures the ratio of an infinitesimal circle around the axis to its diameter to be n. 
The C energy E is then defined as [1, & Gl] 

E:= i(l-r 2 VVV Q r). (2.3) 

o 

This quantity is assumed to represent the line energy density enclosed inside the cylinder. A cylinder, which is a 
two-surface given by u =const and v =const, is said to be trapped, marginally trapped and untrapped if W is 
timelike, null and spacelike, respectively. In terms of the C energy, a cylinder is trapped, marginally trapped and 
untrapped if E > 1/8, E = 1/8 and E < 1/8, respectively. A cylindrical trapping horizon is a hypersurface foliated 
by marginal cylinders. 

The Einstein equations = SnT^ reduce to 

r.uu + %rip 2 u - 27,„r u = 8irrT uu , (2.4) 
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r tVV + 2rip v - 2j^ v r^ v = StttT vv , (2.5) 

r.uv = &TtrT uv , (2.6) 

-2r 2 e- 2 ^(7, tt ,+V,n^) = 8ttT w , (2.7) 

e ^- 2 y(2ril) tUV + r U V>,„ + r^u - r% uv - r, uv - rip )U ip tV ) = 4:irrT zz , (2.8) 

where T^ v is the stress-energy tensor for the matter fields. 



B. Einstein- Rosen waves 

Here we consider vacuum spacetimes in whole- cylinder symmetry. Solutions of the field equations in this case are 
called Einstein- Rosen waves For such a cylindrically symmetric vacuum spacetime, in which the line element is 
given by Eq. (|2.2p . Eq. (|2.6|) reduces to 

r, w = 0, (2.9) 

implying 

r = f(u)+g(v), (2.10) 

where / and g are arbitrary functions. 

If W is spacelike, we can choose / and g by rescaling u and v such that 

(2.11) 



V2 ■ 

Introducing the time and radial coordinates 

v + u , V — u ,„ . 

t= ^jr x= ^r ( } 

we obtain the metric in the form 

ds 2 = e a( 7 -*)(_ dt 2 + dT 2 } + e ~W x 2 dip 2 + e 24> dz 2. ( 2 .13) 

We note that this line element, subject to the Einstein equations below, corresponds to the original Einstein-Rosen 
waves. That is, the Einstein-Rosen paper [l[ deals exclusively with the case where W is spacelike. For convenience, 
we will use the term to refer to any vacuum whole-cylinder symmetric solution of the Einstein equations. 

The nontrivial components of the Einstein equations become the following simple set of partial differential equations: 

-^,tt+V\** + -^,x = 0, (2.14) 
x 

7, x = x(i>% + ^ t ), (2.15) 

7 >t = 2xtP^ jt . (2.16) 

Equation (|2.7p is automatically satisfied due to Eqs. (|2 . 1 5[) and (|2.16[) . The regular axis condition reduces to 7 — > 
as x — ► in this coordinate system. 

If W is timclikc, we can choose / and g by rescaling u and v such that 

r = ^ = , <™> 

Then, we get the metric in the form 

ds 2 = e 2(7-V0(_^2 + ^ 2) + e -2^ t 2 d(p 2 + (2 lg) 

The nontrivial components of the Einstein equations become the following simple set of ordinary differential equations: 

-ip,tt + ip,xx - ji/j,t = 0, (2.19) 

(2.20) 

7, x = 2#,tV',x. (2.21) 
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We can see that the equations are the same as those for spacelike V°r if we exchange t and x. 
If W is null, we can choose / and g such that 

r = u. (2.22) 

The Einstein equations reduce to 

ff> tWI = 0, (2.23) 
1>, v = 0, (2.24) 
%u = u^ 2 u . (2.25) 

ip = i>{u) follows from the above. 



C. Self-similar spacetimes in whole-cylinder symmetry 

We now consider the case where the spacetime is self-similar as well as cylindrically symmetric: there is no a priori 
reason to suppose that this will lead to trivial solutions only. In other words, we assume that the spacetime admits a 
vector field v which satisfies the following equation 

£ v flV = 2 9nu, (2.26) 

where £ v denotes the Lie derivative along v. We refer to the vector field v and Eq. (|2.26p as the homothetic vector 
and the homothetic equation, respectively. 
We assume that v has the following form: 

Q Q 

v = a(u,v)—+P(u,v)—, (2.27) 
ou ov 

where v is then assumed to be orthogonal to the cylinders of symmetry. We will refer to this as a cylindrical homothetic 
vector. We will mention the limitation of this assumption later. Then, the homothetic equations (|2.26[) yield a = a(u) 
and (3 = f3(v). We can then gencrically make the coordinate transformation u = u(u) and v = v(v) satisfying the 
following relations: 

du dv 

a(u)^=2u and (3{v)— = 2v. (2.28) 
du dv 

Then, we can have 

v = 2u— + 2v—, (2.29) 
ou ov 

where here and hereafter we omit bars for simplicity. The homothetic vector is timclikc, spacelikc and null if uv is 
positive, negative and zero, respectively. If a = or (3 = 0, the homothetic vector is null and we do not consider these 
special cases here. 

It is then straightforward to show that the homothetic equations (|2 . 26[) lead to 

e 2 ^ = Me 2p ("\ r=\u\S(v), e 2 ~< = e 2G ^ (2.30) 

where 77 = v/u and P, S and G are arbitrary functions. 

Thus, we obtain the following standard form of the metric for whole-cylindrically symmetric self-similar spacetimes: 

ds 2 - -2e 2 °W-* p W\ u \-i dudv + e -mn) \ U \S 2 ( V ) dip 2 +e 2P ^\u\dz 2 . (2.31) 



We can substitute the above form into the Einstein equations (|2.4|) - (|2.8|) and obtain a set of ordinary differential 
equations for P, G and S. Fortunately, for the vacuum case, we can greatly simplify the system. 
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III. EINSTEIN-ROSEN WAVES WITH A CYLINDRICAL HOMOTHETIC VECTOR 



A. One-parameter family of solutions 



We consider self-similar vacuum solutions in this section. It is not trivial that the choice of / and g adopted 
in Sec. HIB is compatible with the self-similarity introduced in Sec. QTJI In fact, for r to be compatible with the 
self-similarity, i.e., Eqs. (|2.30[) . / and g introduced in Eq. (|2. 10[) for the vacuum solution must satisfy 



/ = CiU, and g = C%v, 



(3.1) 



where C\ and C2 are arbitrary constants. (We note that a trivial addition of a constant may be required to obtain the 
above form.) Hence, the choices of / and g given by Eqs. (|2.1ip and (|2.17|) are both compatible with self-similarity. 

For the moment we restrict ourselves to the case where r has a spacelikc gradient. Now we can adopt (t, x) 
coordinates and by combining Eqs. (|2. 13|) . (|2.30p and (|2.3ip . we find the metric in the following form: 



where 



(We recall that x 
equations. Noting 



ds 2 = e 2r«)-2*«) a .-l ( _ £K 2 + ^2) + e -2»(Q xd/p 2 + ^{6) xdz ^ 



(3.2) 



(3.3) 



> 0.) In this case, the Einstein equations reduce to a set of simple ordinary differential 



1 



¥(0 + 2 In a:, 



j, = 

7 = r(0, 

Eqs. (|2.14|) . (|2. 15[) and (|2.16|) reduce to the following ordinary differential equations: 

(£ 2 - !)¥" + £¥' = 0, 



r' = -2 ( £¥' - 



-¥' 
¥', 



where the prime denotes the derivative with respect to £. From Eqs. (|3.7|) and (|3.8p . we obtain 



1 



Therefore, ^ 2 > 1 and 



¥' = ±- 



1 



This satisfies Eq. (|3.6p . We can integrate the above equation and obtain 

1 



¥ = - In 
2 



¥r 



where ¥o is a constant of integration. T is obtained by substituting Eq. (|3.10p into Eq. 
resultant equation. The result is 



r = Im 

2 



± 1 



+ A, 



(3.4) 
(3.5) 

(3.6) 
(3.7) 

(3.8) 
(3.9) 
(3.10) 

(3.11) 

and integrating the 
(3.12) 



where A is a constant of integration. Getting back to the original metric functions tp and 7, the solution is given by 

V; = iln|e± V / F _ T| + ¥o + 7 ln ^ (3-13) 



In 



± 1 



(3.14) 
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We can assume t > because the flip of the sign of t corresponds to the other branch of solutions. Note that we can 
set = by absorbing it into the coordinates as follows: 

t = e - 2 *°t, £ = xe" 2 * , (p = ip,z = e 2 *°z. (3.15) 

Up to this gauge parameter, the solutions are parametrized by A. 

For the timelike W case, we can obtain the solution just by exchanging t and x in the solution for the spacelikc 
W case. This corresponds to the replacement of £ with The solution is therefore given by 

i> = \ In If" 1 ± VZ- 2 - 1| + *o + \ In |*|, (3.16) 
1 

7 = 2 

where the solution is valid for f 2 < 1. We can assume x > in this case. 



± 1 



+ A, 



(3.17) 



B. Flatness and topology of the solutions 

As indicated in the introduction, our aim in deriving the solutions of the previous subsection is to study the 
simplest case of self-similar cylindrical collapse. However it transpires that these solutions do not represent collapsing 
gravitational waves. The solutions are flat everywhere except along the axis, and thus the solutions correspond 
either to part of Minkowski spacetime, or to a line conical singularity in flat spacetime: the assumption of self- 
similarity of the gravitational waves rules out any other possibility. This is a nontrivial result regarding self-similar 
cylindrical collapse. Furthermore, the form of the flat spacetime metric that emerges demonstrates explicitly that 
there are cylindrical trapping horizons in Minkowski spacetime, and that C— energy is (i) nonunique and (ii) nonzero 
in Minkowski spacetime. This seriously undermines the interpretation of C— energy as the gravitational energy of a 
cylindrical spacetime. 

We can explicitly show that all coordinate components of the Riemann curvature tensor vanish for the obtained 
solutions. This means that the solutions are flat. Here we show that these solutions are indeed part of the Minkowski 
spacetimcs. 



Untrapped case 



We here assume that W is spacelike, i.e., (i, x) corresponds to an untrapped cylinder. The regular axis condition 
implies r — > as f — > oo. From Eq. (|3 . 1 2|) . this is possible only for the upper-sign solution with A = 0. Hence, the 
solution becomes 



1 



7 



In 



1 



2 VV^T 



(3.18) 



The upper-sign solution for the different choice of A gives a conical singularity with the ratio of an infinitesimal circle's 
circumference to its diameter 7re _A rather than tt. We first concentrate on the upper-sign solution with A = below. 
Choosing the coordinates where vPq = 0, we can write the metric as 



ds 2 



1 



,(-dt 2 + dx 2 ) + \t T Vt 2 - x 2 \d<p 2 + \t±y/t 2 -x 2 \dz 2 



where the solution is valid only for t 2 > x 2 . Assuming t > 0, through the coordinate transformations 



T 2 = t + \/t 2 - x 2 and X 2 = t~ \f¥ 



(3.19) 



(3.20) 



t 



T 2 +X< 



and x = TX, 



we obtain the following metric for the upper-sign: 

ds 2 = -dT 2 + dX 2 + X 2 dip 2 + T 2 dz 2 , 



(3.21) 
(3.22) 
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where < X < T by construction. Through another coordinate transformation 

t = T cosh z, C = T sinh z, p = X cos tp and q = X sin (3.23) 

we finally obtain the usual Minkowski spacetime in the standard Cartesian coordinates 

ds 2 = -dr 2 + dp 2 + dq 2 + d( 2 , (3.24) 

where r 2 > p 2 + q 2 + £ 2 and hence the solution covers the inside of the light cone r 2 = p 2 + q 2 + £ 2 . This region 
is shown as the dark shaded disk in Fig. [T] where the constant r spacelike hypersurface is plotted. The z-axis, i.e. 
x = in (t,x,ip,z) coordinates, is transformed to the £-axis, i.e., p = q = in (r,p,q,£) coordinates. The null 
hypersurface t 2 = x 2 in (t, x, ip, z) coordinates is transformed to the light cone r 2 = p 2 + q 2 + £ 2 in the standard 
Cartesian coordinates. 




FIG. 1: The constant r spacelike hypersurface is shown in the standard Cartesian coordinates (r,p,q,£). The dark shaded 
region is untrapped, while the light shaded region is trapped. The circle shows the light cone, which is a cylindrical trapping 
horizon. The region which is unshaded is not described by the cylindrical vacuum flat solutions. The dashed lines denote the 
timelike planes which are identified with ( = 0. 



For the lower-sign solution, the metric is written as 

ds 2 = -dT 2 + dX 2 + T 2 df 2 + X 2 dz 2 , (3.25) 
where < X < T. Through another coordinate transformation 

t = T cosh ip, £ = T sinh (p, p = X cos z and q = X sin z, (3.26) 

we finally obtain the metric 

ds 2 = -dr 2 + dp 2 + dq 2 + dC, 2 , (3.27) 

where r 2 > p 2 + q 2 + Q 2 . This is also the Minkowski spacetime in the standard Cartesian coordinates but with 
nontrivial topology. In the original (t, x, ip, z) coordinates, (p = 2ir is identified with <p = 0. This results in the 
identification between the two timclike hypcrsurfaccs £ = and £ = Vt, where V = tanh27r. The latter timclike 
hypersurface is shown as a dashed line in Fig. [TJ More precisely, the point (r, p, q, 0) is identified with the point 
(t/VI - V 2 ,p, q, tV/VI - V 2 ). Then, there appears a timelike geodesic with infinite spatial acceleration in an ap- 
proach to the spacclikc line r — £ = 0, which will be described in detail in the appendix. The appendix also serves to 
clarify the nature of the topological identifications made here. On the other hand, a circle identification has not been 
imposed on z, implying that the two-dimensional pq plane consists of covering planes folded infinitely many times. 
The z-axis is again transformed to the C-axis, i.e., p = q = 0. 

We note that here the (original) cylinders of symmetry (t, x) = (to, xq) with to, xo both constant have the following 
representation in Minkowski coordinates (r,p, q, £): 



t 2 - C 2 = *o + V^o-^O' p 2 +q 2 = t - \jtl~x 2 



(3.28) 
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2. Trapped case 

We then assume that W is timelike, i.e., (t, x) corresponds to a trapped cylinder. Also for this case, choosing 
^0 = and A = 0, we obtain 

ds 2 = (-dt 2 + dx 2 ) + \xt \/x 2 -t 2 \d(p 2 + \x± ^x 2 -t 2 \dz 2 , (3.29) 

2V x 2 — t 2 

where the solution is valid only for x 2 > t 2 . By implementing coordinate transformations similar to those used in the 
untrapped case, we obtain for the lower-sign solution 

ds 2 = -dT 2 + dX 2 + X 2 dip 2 + T 2 dz 2 , (3.30) 

where < T < X. This is identical with the upper-sign solution for the untrapped region and hence transformed 
to the standard Cartesian coordinates (r,p, q, (), where the solution covers the region £ 2 < t 2 < p 2 + q 2 + ( 2 . This 
corresponds to the intersection of the outside of light cone t 2 = p 2 + q 2 + ( 2 and the timelike portion sandwiched by 
two planes r = ±£. This region is shown as a light shaded region in Fig. [1] The topology of the spacetime is trivial. 
For the upper-sign solution, we obtain 

ds 2 = -dT 2 +dX 2 +T 2 dip 2 + X 2 dz 2 , (3.31) 

where < T < X. This is identical with the lower-sign solution for the untrapped region and hence transformed to 
the standard Cartesian coordinates (r,p, q, £), where the solution covers the region ( 2 < t 2 < p 2 + q 2 + ( 2 7 i.e., the 
intersection of the outside of the light cone and the timclike portion sandwiched by two planes r = ±£. We should 
note that the topology is nontrivial because (p is circularly identified while z is not. 

We note that in the trapped case the (original) cylinders of symmetry (t,x) = (to,Xo) with to,xo both constant 
have the following representation in Minkowski coordinates (r,p,q,Q: 

t 2 - ( 2 = xq - yx 2 -* 2 , p 2 +q 2 = x Q + ^x 2 -t 2 . (3.32) 



Hence, for A = 0, the union of the upper-sign solution for the untrapped region and the lower-sign solution for the 
trapped region describes the timelike portion of the Minkowski spacetime sandwiched by two light planes when the 
two solutions are matched on the light cone r 2 = p 2 + q 2 + ( 2 . 

The union of the lower-sign solutions for the untrapped region and the upper-sign solution for the trapped region 
also describes the same region of the Minkowski spacetime when they are matched on the light cone r 2 = p 2 + q 2 + £ 2 
but with nontrivial topology. The solutions with nonvanishing A will have an additional conical singularity. 

This solution without the conical singularity is quite analogous to the Milne universe solution, which is also part 
of the Minkowski spacetime. However, the present solution is somewhat different from the Milne universe in the 
following respect. Recall that observers with constant spatial coordinates run radially outward with a constant speed 
and they do so homogeneously and isotropically in the Minkowski spacetime. On the other hand, in the present 
solution, observers with constant spatial coordinates run with a constant speed but only in the direction parallel to 
the axis and do not in the two perpendicular ones. We shall call the present solution the cylindrical Milne solution in 
this paper. 

C. C energy and trapping horizon in the Minkowski spacetime 

C— energy was introduced in Q as a tool with which cylindrical spacetimes may be analyzed. It has several 
interesting and useful features: It is covariant and is associated with a conserved flux vector; it has the correct 
Newtonian limit, the mass per specific length of the cylinder [To| ; it is propagated by Einstein- Rosen waves. Thus 
it is a candidate for "the energy of whole- cylinder- symmetric spacetimes^ (the phrase appears in quotation marks 
in 0], p. 251) and a later study refers to C— energy as "gravitational energy per specific length" As a particular 
application, C— energy has been used to investigate the fate of an infinitesimally thin cylindrical shell composed 
of counter-rotating dust particles by Apostolatos and Thorne 0], and later by one of the present author (KN) and 
his collaborators [Tj]: it should be stressed that the conclusions in these two papers do not agree with each other 
although both rely on the properties of the C— energy. A further criterion that should be satisfied by a candidate for 
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gravitational energy of any form is that it should vanish in the absence of a gravitational field, i.e. in flat spacctime. It 
transpires however that the cylindrical representations of flat spacetime we have found above show that the C— energy 
does not always vanish in this case. 

For the Einstein-Rosen waves written in the forms of Eqs. (|2.13[) and (|2.18p , the C energy reduces to the following 
simple forms: 

E=\{l-e- 2 ^) and E=\{l + e- 2 ^), (3.33) 
o o 

respectively. For simplicity, we discuss the cylindrical Milne solution with trivial topology, which is given by pasting 
the upper-sign solution of Eq. (|3 . 19[) and the lower-sign solution of Eq. (|3.29p on the null hypersurface t 2 = x 2 . For 
these metrics, we obtain respectively 



E=lt -vw^ and ix + VS^f 

8 t + s/i 2 - x 2 8 x - Vz 2 - t 2 



This is rewritten in both cases as 

1 p 2 + q 2 

E = g r 2 - C 2 ' ( ] 

in terms of the standard Cartesian coordinates. Since these metrics arc those for part of the Minkowski spacetime, it 
can have nonvanishing C energy. 

However, if we write the metric of the Minkowski spacetime in the standard cylindrical coordinates, we have 
7 = ip = in Eq. (|2.13[) and hence E = 0. This result questions the physical interpretation of the C energy. Indeed, 
the trick is in the choice of the two commuting Killing vectors, or equivalently, the choice of the cylinders. If we take 
(£(i),£(2)) = (d/dip,d/dz), then p, £ and r are calculated as 



P 



\t-Vt 2 -x 2 \=p 2 + q 2 , (3.36) 



( 2 = \t+y/t 2 -x 2 \=r 2 -C 2 , (3.37) 
r 2 = x 2 = (r 2 -Q 2 ){p 2 + q 2 ), (3.38) 

where t 2 > x 2 or r 2 > p 2 + q 2 + ( 2 . Hence, the regular axis condition is satisfied and we obtain nontrivial C energy. 
Instead, if we take (£(i), £(2)) = (9/ dip, d/dQ, then we have p 2 = p 2 + q 2 , £ = 1 and r 2 = p 2 + q 2 and obtain vanishing 
C energy. We should also note that the z-axis in the former is transformed to £-axis in the latter. This means that 
the definition of C energy is ambiguous for the same axis in the same cylindrically symmetric spacetime unless a pair 
of two commuting Killing vectors are fully specified. Thorne [3] has noted the lack of uniqueness in the definition of C 
energy in the case of unpolarized cylindrical spacetimes, for which the Killing vectors (£(i),£(2)) are not orthogonal. 
Nonuniquencss in the unpolarized case is related to the loss of invariance of the spacetime under reflections through 
any plane cither containing the axis or perpendicular to it. However as we see in the present case, nonuniqueness can 
remain even in the polarized case when there is more than one choice of the azimuthal and translational Killing fields. 

This clearly gives rise to a question about the interpretation of C energy as 11 gravitational energy per specific 
length" [13], given that it may be nonzero in the absence of a gravitational field. But as we have pointed out, C 
energy has many attractive and useful features, and so perhaps the most natural question to ask at this point is if 
there exists an alternative definition that would have the additional feature of vanishing for any cylindrical slicing of 
Minkowski spacetime. We hope to address this question in future work. 

Although the uniqueness of the C energy may be recovered by specifying the pair of Killing vectors, it is still true 
that the null hypersurface t 2 = x 2 in the original coordinates or the light cone r 2 = p 2 + q 2 + C 2 in the standard 
Cartesian coordinates gives a cylindrical trapping horizon. The inside of the light cone r 2 > p 2 + q 2 + ( 2 is untrapped, 
while the outside of the light cone, i.e., £ 2 < t 2 < p 2 + q 2 + ( 2 is trapped. The constant r hypersurfaces, given by 
Eq. (|3.38[) . are shown in Fig. [2] on the constant r hypersurface. On the other hand, with the pair (d/d(p,d/d() of 
Killing vectors, the constant r hypersurfaces given by 

p 2 + q 2 = r 2 (3.39) 

are all timelike. See also [13] for trapped surfaces in the Minkowski spacetime. 

In cylindrically symmetric spacetimes, a trapping horizon is defined as a hypersurface foliated by marginally trapped 
cylinders and hence will not be closed in general. Since they are not closed, it will not imply the existence of spacetime 
singularity. Thus, this example is a lesson that we cannot reasonably identify a trapping horizon with a black hole 
horizon for cylindrically symmetric spacetimes [ToL [l2| . On the other hand, we could have closed trapped surfaces if 
we change the identification. In that case, we may encounter a sort of singularity as described in the appendix. 
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p,q 

FIG. 2: The constant r spacelike hypersurface is shown in the standard Cartesian coordinates (r,p,q,(), where r is chosen 
to be unity. The short-dashed, dashed and long dashed lines denote the intersection of the constant r hypersurfaces given by 
Eq. (|3.38p with the constant r hypersurface. 



IV. MORE GENERAL SELF-SIMILAR EINSTEIN-ROSEN WAVES 

In the previous two sections, we showed that the only self-similar cylindrically symmetric vacuum spacetimes 
comprise flat spacetimes, possibly with line conical singularities along the axis. These are trivial examples of self-similar 
spacetimes, and so while these may act as (intermediate) asymptotic endstates of certain more general cylindrical 
configurations - and indeed a line conical singularity in an otherwise flat spacetime can result from the complete 
collapse of cylindrical null dust - it is clear that this does not lend any weight to the self-similarity hypothesis in 
cylindrical symmetry It transpires however that part of the reason for this is that we have looked only at a quite 
restrictive class of self-similar cylindrical spacetimes. Dropping the assumption that the homothetic vector field is 
orthogonal to the cylinders of symmetry yields some interesting results which we describe here. 

A. Two-parameter family of solutions 

In this section, we concentrate on the spacelike W case, where the region is untrapped. In this case, Eqs. (|2.14l) . 
(I2.15|) and (|2.16[) give a complete set of governing equations. Among them, Eq. (|2.14[) gives the dynamics of ijj and 
the other equations determine the derivatives of 7. For self-similar solutions with a cylindrical homothetic vector, tp 
is given by Eq. (|3.4|) . Note that each of two terms on the right-hand side of Eq. (|3.4j) is a solution of Eq. (|2.14|) - 
provided that the first term 'J is taken to be a solution of Eq. (|3.6p . Since Eq. (|2 . 14[) is linear, this means that the 
arbitrary linear combination of the two terms gives a solution of Eq. (|2.14p . Moreover, we can assume the similar 
form also for 7 from Eq. (|3.5[) . From this consideration, we assume the following form for -0 and 7: 

tJj = a*(Q + b±\n\x\, (4.1) 
7 = cf(0+^m|4 (4.2) 

From Eqs. (|2.14p . (|2.15p and (|2.16p . we get the following ordinary differential equations for \P and T: 

(£ 2 _ + ' = 0, (4.3) 
+ 2 + a 2 *' 2 , (4.4) 

cf ' = -2 ( a£*' - b -) of. (4.5) 
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From Eq. (|4.3|) , we get 



(4.6) 



We can always assume ^1 = 1/2 because of the factor a in Eq. (|4.ip . Eliminating f from Eqs. ()4.4p and (|4.5|) . using 
Eq. (|4.6|) and putting = 1/2, we obtain 



sign(l -£ 2 )a 2 + b 2 = 2d. 



Then, for £ 2 > 1, the solution is given by 



1 



V- = a- In |£ + VC 2 - 1| +6-lna; + *o, 



1 



a— 
2 



- a l nv /e 2 -l + 61n|e+ v^ 2 ~ 1| 



d-ln|x| 



(4.7) 

(4.8) 
(4.9) 



where we omit the lower-sign solution because the sign can be absorbed into the sign of a. 

To be more specific, we assume that the axis is regular or conically singular at least, which implies that 7 approaches 
a finite value for t > and x — > 0. This condition strongly restricts the parameters. This limit corresponds to £ — > 00 
and x — > 0, where £ _1 and a; approach zero independently. Hence, the condition on the axis implies a = b, d = and 
a finite value for Tq- Putting a = b = 2k and Tq = A, we obtain the following solution with a regular or conically 
singular axis: 



2k z In 



ln(£ + V^ 2 - 1) + In M 



1 



+ A, 



(4.10) 
(4.11) 



where V&o is eliminated in use of the scaling freedom of t and x. These solutions are parametrized by k and A. Note 
that the solution is self-similar in the sense discussed so far if k = 1/2 and the axis is regular if and only if A = 0. As 
we will see later, the spacetime is nonflat except for k = and 1/2. 

For later convenience, we write down the line element explicitly both in (t, x, ip, z) and (T, X, ip, z) coordinates as, 



If 1 A /+2 _ „2\2k(2k-1) 

= * + L„.., e 2 \~dt 2 +dx 2 ) + 



2 i^(t 2 -x 2 ) 



2\2k 2 



(t + Vt 2 - X 2 f> 



-dip 2 + (t + Vt 2 - x 2 ) 2K dz 2 



ji4k;(2k;-1) x 2 

v_ YW -i e 2X (~ dT2 + dX ") + ^-ry dp 2 + T^dz 2 , 



(T 2 - X< 



(4.12) 



where T and X are given by Eq. (|3.20[) . The original domain of the solution is given by0<x<£<oo and this is 
mapped to < X < T < 00. It is clear that for k = the solution reduces to the Minkowski spacetime with a regular 
(A = 0) or a conically singular (A 7^ 0) axis. Moreover, we can easily find that for k = —1/2 this reduces to a Kasner 
solution with a regular (A = 0) or a conically singular (A 7^ 0) axis. 



B. Non cylindrical homothetic vector 



To make it clear whether the solutions obtained above have some kind of self-similarity, we shall consider the scaling 
transformation i = At, x = Ax, <p = p and z = z. Through this transformation, tjj and 7 transform as follows 



ijj = ijj + 2k In A and 7 = 7. 
Then, the metric components g tt , g xx , g vv and g zz transform as follows: 

g tt = A 2(1 ~ K) 9tu 9xx = A 2il - K) g xx , g lpv =A 2 ^-^g tptp and g„ = A 2K g zz . 
Therefore, for k ^ 1, if we define a vector field v as 

1 8 1 d 

1 — k Ot 1 — k dx 1 



(4.13) 
(4.14) 

(4.15) 
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we can write 

2re 

C v g t t = 2g u , C v g xx =2g xx , C v g vtp = 2g tptp and C v g zz = g zz . (4.16) 

I — K 

It is clear that if and only if re = 1/2, the vector field v is a homothetic vector. 

However, if we instead scale the coordinates as t = At, x = Ax, ip = <p and z = A x ~ 2k z, the metric components 
transform as follows: 

fa = A^-^g^. (4.17) 
For re 7^ 1, this implies that for the vector field w given by 

w — 1 t 9 + 1 x 9 + 1 ~ 2K z 9 (4 18) 
1 — k dt 1 — k dx 1 — re dz 

we obtain 

£wflV = ^9nv (4.19) 

Therefore, w is a homothetic vector and hence the spacetime described by this solution is self-similar. The homothetic 
vector w is not cylindrical if re ^ 1/2. 
For re = 1, if we instead define W as 

d d d 

W:=t^-+x^-z^-, (4.20) 
at ox oz 

we obtain 

Cwgab = 0, (4.21) 
and therefore W is a Killing vector. This is independent from the azimuthal and the translational Killing vectors. 



V. PHYSICAL INTERPRETATION OF SELF-SIMILAR EINSTEIN-ROSEN WAVES 

In the previous section we have derived a two-parameter family of self-similar solutions, where the homothetic vector 
is not cylindrical in general. As we have seen, the two-parameter family of solutions includes a Minkowski (k = 0), 
cylindrical Milne (k = 1/2) and Kasner (k = —1/2) solutions. In this section we present the physical interpretation of 
the family of solutions for general values of re, based on the analysis of spacetime singularities, infinities and analytical 
extensions. We show that these solutions describe physically interesting nonlinear phenomena of gravitational waves. 



A. Curvature invariant, C— energy and areal radius 



For the two-parameter family of solutions obtained here, the Kretschmann invariant / := R^P" R^ vpa 
to give 



is calculated 



I = 2 4+8K re 2 (l - 2re) 2 e" 4A (t + ^/W^)- 2 ^ ~ 2K+1 \i 2 - x 2 ) 4 " 3/2 
(1 + 2re)(l - K)t + (2 - re+ 2 K 2 )\/t 2 - a 
2 6 re 2 (l - 2re) 2 e- 4A T- 4 ( 4K2 - 2K+1 )(T 2 - x 2 ) 8k2 ~ 3 [3T 2 - (4re 2 - 2re + 1)X 2 



(5.1) 



This is identically zero and hence the spacetime is flat if and only if re = or 1/2. For T 2 = X 2 the invariant is 
diverging if < re 2 < 3/8 and re 2 ^ 1/4, while it is finite if re 2 > 3/8 or re = 0, ±1/2. The invariant is finite at the 
z-axis orx = 0if0<t<oo. The invariant at the axis is vanishing even for 0<t<ooifre = 0, re=l/2orre>l, 
while it is diverging at£ = 0ifre<0,0<re< 1/2 and 1/2 < re < 1. Only for re = 1, it is nonzero and finite for 
< t < oo on the axis. Since the solution is vacuum, the Riemann tensor reduces to the Weyl tensor. Hence, in a 
intuitive sense, we can regard the Kretschmann invariant as the field strength of the pure gravitational field and/or 
gravitational waves, although this may be negative for re < —1/2. 
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The C energy for this solution is calculated to give 



E 



1 - e 



-2,\ 



2 \^/t 2 -x 2 





i 




r; 


~~ 8 



-2A 



T 2 -X s 
2^ 



(5.2) 



£ = (1 



-2A 



)/8 at the axis x 



0. This is nonzero if and only if A ^ 0, i.e., at the conical singularity. Whether 
A is zero or not, E = 1/8 for t 2 = x 2 if only k / 0, suggesting a cylindrical trapping horizon. For n = 0, the C 
energy is identically given by E = (1 — e~ 2A )/8 and hence constant in the whole region described by this solution. As 
already mentioned, E is vanishing for k = A = 0, but is not for k = 1/2 and A = 0, although both correspond to a 
flat geometry. The C energy is unchanged by flipping the sign of k. 
The areal radius r of the cylinder can be calculated as 



r 2 = x 2 = T 2 X 2 . 



(5.3) 



Hence, the r =const surface is given by a hyperboloid in the TX plane. 



B. Analytical extension and global structure 

We recall that the original domain of the solutions is mapped to < X < T < oo in TX plane. The "event" 
t = x = 0, or equivalently, T = X = seems to be a s.p. (scalar polynomial) curvature singularity because of the 
divergence of the curvature invariant if K < 1 and n =/= 0, 1/2. But this is a subtle issue. The proper time a along a 
curve on which X vanishes and z is constant is given by 

a = e A J T~ 2K+1 dT. (5.4) 

We can easily see that a is finite in an approach to T = if and only if k < 1. If and only if k > 1, then T = X = 
is at a timelike infinity and hence not within the physical spacetime. On the other hand, (T, X) = (oo, 0) is at a 
timelike infinity for n < 1, while it is a spacetime singularity for k > 1. 

For t 2 — x 2 , or equivalently, T 2 = X 2 , the curvature invariant is diverging as long as < k 2 < 3/8 and k 2 ^ 1/4. 
We introduce the following null coordinates: 

U = T-X and V = T + X. (5.5) 

Then, we have 

,2 2X (V + U\ iK[2K - 1] dUdV f V + U\ 2{1 - 2K) fV-U\ 2 J a fV + U\ 4 * 2 

The original domain of the solutions is mapped to0<£/<y<oo. Here we consider radial null geodesies along 
which U , ip and z are constant. The geodesic equation for these null geodesies is given by 



d 

du> 



(V + (/) 4k ( 2k - 1 ) dV 



(VU)^ 2 - 1 da; 



0, (5.7) 



where U =const and u> is an affine parameter. Thus, we can find that V = oo is a null infinity for any n. We can 
similarly find that U = can be reached in a finite affine length along null geodesies with ^=const for < k 2 < 1/2, 
while it is a null infinity for k 2 > 1/2. Below we discuss the cases k 2 = 1/2,0<k 2 <1/2 and k 2 > 1/2, separately. 



1. k 2 = 1/2 

First, we consider the case of k 2 = 1/2, and introduce 

U = e u and V = e v . (5.8) 

Then, the metric (15.61) becomes 



2(2=fV2) / v I e u\ 2(l=F\/2) /„u_„u\ 2 / e v i e u \ ±2^2 

dg > = , j dudy + ^ 2 + dz 2 (5>g) 
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for« = ±l/V2. The original domain is mapped to— oo < u < v < oo and this describes a maximal extension because 
v = oo and u = — oo both correspond to null infinities. In the above coordinate system, there seems to be no singular 
point. But as already shown, the "event" u = v = — oo or T = X = is a spacetime singularity. The conformal 
diagram of the solution is given in Fig. [3(d)] There is no cylindrical trapping horizon and the whole spacetime is 
untrapped. 



2. < k 2 < 1/2 

Next, we consider the case of < k 2 < 1/2. In this case, we introduce u and v as 

U = u n and V = v n , (5.10) 
where n := 1/[2(1 - 2k 2 )] > 0. Then, we have 



*• . -pa-*™--" (^f- 2 " (K_£) 2 v + (±^) 4 V. (5 .n) 

The original domain is mapped to < u < v < oo. As we have shown, v — oo is a null infinity, while u = is finite. 
For < k 2 < 3/8 and k 2 ^ 1/4, the Kretschmann invariant diverges at u = and hence tt = corresponds to a 
null singularity. For 3/8 < k 2 < 1/2, u = is an at least C 2 extendible null hypersurface, which coincides with a 
cylindrical trapping horizon, and we can discuss the extension beyond this hypersurface. To examine the affine length 
of the radial null geodesic with u = 0, we should consider 



d 

du> 



(V + mMto-i)* 
dw 



= 0, (5.12) 



instead of Eq. (|5.7[) . We can then find that the affine length is infinite to v = oo and finite to v = even along the 
null geodesic with u = for 3/8 < n 2 < 1/2. 

If and only if n is a natural number, the extension beyond this surface can be analytic and we can naturally extend 
the spacetime by Eq. (I5.10|) . The following discussion depends on whether n is odd or even. 

If n = 21 + 1 (/ = 1, 2, ■ ■ ■), the maximally extended domain is given by {0 < v < oo and — v < u < v}. On 
the surface v = —u, we have T — and hence the Kretschmann invariant diverges. So this surface corresponds to a 
spacelike singularity. There is a cylindrical trapping horizon on u = 0. The region {0 < u < v < oo} is untrapped, 
while the region {0 < — u < v < oo} is trapped. 

If n = 21 (I = 1, 2, • • •), the maximally extended spacetime is apparently given by — oo < u < v < oo. u = — oo is a 
null infinity, v = — u = oo is a spacclikc infinity, u = v = is a spacetime singularity, v = u = ±oo are both timclike 
infinities. It is interesting to see the surface v = —u. Noting T = [v n + u n )/2 and X = (v n — u n )/2, we find that 
the Kretschmann invariant is finite there except for u = v = 0, while the areal radius r vanishes. It turns out that 
we need to pay close attention to this surface. To get an insight into this surface, we introduce t and x coordinates, 
where u = t — x and v = t + x, so that we should focus on the surface t = 0. Near this surface, the metric can be 
written as 

ds 2 ~ [2(1 - 2n 2 )]~ 2 x 2 < 2K -^/^- 2 ^ [e 2X (-dt 2 + dx 2 ) + t 2 d^ 2 } + x^^-^Uz 2 . 

It follows from the identification between ip — and tp = 2ir that there is a timclike geodesic on the ttp plane with 
infinite spatial acceleration, as shown in the appendix, and hence the spacetime is geodesically incomplete in a sense 
that there is a geodesic which cannot be uniquely extended. Thus, there is no analytical extension beyond this surface. 
Then, the structure of the resultant spacetime is similar to the n = 21 + 1 case, except for that "singularity" on the 
surface t = is only "quasiregular" in the sense of Ellis and Schmidt [3(| ■ 

It should be noted that it is impossible to analytically extend the spacetime even beyond u = even for 3/8 < 
k 2 < 1/2 if n — 1/ [2(1 — 2k 2 )] is not an integer. For this case, the functions U and V are at least twice diffcrcntiablc 
but not C°° with respect to u at u — 0. The spacetime admits at most extension beyond the null surface u = 0, 
where [n] is the largest integer which is no greater than n. 



The conformal diagrams of the solutions for different values of n are given in Figs. 3(a) and 3(b) As we can see 
in these figures, the case of integer n is particularly intriguing in the context of gravitational collapse because these 
solutions are self-similar, describe the collapse of gravitational waves and admit nonsingular initial data on a spacelike 
Cauchy surface containing both trapped and untrapped cylinders. 
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3. k > 1/2 

Also in this case, we introduce u and v as 

U = (-u) n and V=(-v) n , (5.13) 

where n = 1/ [2(1 — 2k 2 )] < and we put the negative signs to keep the increase of u and v corresponding to that of 
U and V, respectively. Then, the line element is given by exactly the same form as Eq. (|5.11l) . The original domain 
is mapped to — oo < u < v < 0. Since both u = — oo and v = correspond to null infinities, the original domain 
describes the whole spacetime. Noting T = [(-«)" + (-u) n }/2 and X = [{-v) n - (-u) n ]/2, we find that T = X = 
corresponds to v — u = — oo, while (T, X) = (oo,0) corresponds to v = u = 0. The conformal diagrams of the 
solutions for different values of k arc therefore given in Figs. [3(d)] [3(c)] and [3(f)] There is no cylindrical trapping 
horizon and the whole spacetime is untrapped in this case. 



In summary of this section, the two-parameter family of solutions describe a variety of global structures. They are 
classified in terms of k, or cquivalently, n — [2(1 — 2k 2 )] -1 . For k = 0, 1/2 and —1/2, the solution reduces to the 
Minkowski, the cylindrical Milne and the Kasner solutions, respectively. If < k 2 < 1/4 or 1/4 < k 2 < 3/8, the 
spacetime describes the interior of the exploding (imploding) cylindrical shell of gravitational waves. For 3/8 < k 2 < 
1/2, we have the following three cases: if n = 3, 5, 7, ■ ■ •, the spacetime describes the collapsing gravitational waves 
to a spacelikc singularity or exploding gravitational waves from a spacclikc curvature singularity; if n = 2,4, 6,- • ■, 
the spacetime structure is quite similar to the odd n case but the spacelike curvature singularity is replaced by a 
quasiregular one; if n > 2 is not an integer, the spacetime does not admit an analytic extension beyond the null 
surface. For k 2 > 1/2 and k ^ 1, the conformal diagram is similar to the Minkowski one except for that a singularity 
replaces a timclikc infinity. For k = 1, the conformal diagram is similar to the Minkowski one and the whole spacetime 
is regular. 

It is interesting to note that the present analysis proceeds quite analogously to Wang's [IH for cylindrically sym- 
metric self-similar solutions with a massless scalar field in (3+l)-dimensions and hence Hirschmann's et al. [3l| for 
circularly symmetric self-similar solutions with a massless scalar field in (2+l)-dimensions although the system and 
the result are both different in detail. This follows from the fact that the governing system of partial differential equa- 
tions are quite similar for these systems. It is also interesting to note that the present solution generically involves 
singularities, because one does not usually get a singularity from the collapse of cylindrical waves (global regularity 
of Einstein- Rosen waves: see e.g. Ashtckar ct al. @). In order to generate the singularity, initial data on the Cauchy 
surface must be non asymptotically flat. So we have a sort of converse result: if we allow non asymptotically flat initial 
data, then a singularity can form solely from the collapse of gravitational waves. Furthermore, if trapped surfaces are 
initially absent, then they cannot form (see Thorne Q). In the relevant cases (Fig. |3(b"J[ integer n > 2), we have a 
nonsingular Cauchy surface which already contains trapped cylinders. 



VI. CONCLUSIONS 



In this paper, we have studied self-similar vacuum spacetimes in whole-cylinder symmetry: self-similar Einstein- 
Rosen waves. The primary motivation was to determine possible (intermediate) asymptotic endstates for more general 
Einstein-Rosen waves, and for other cylindrical systems. That is, we wish to study the self-similarity hypothesis in 
cylindrical symmetry, especially in the context of g ravitational collapse. There is a considerable body of evidence 
for the hypothesis in spherical collapse (see e.g. [jj|) and there is also evidence for the hypothesis in the context of 
cosmological models: Hewitt et al. have shown that among a class of cylindrical inhomogeneous cosmological models, 



there are self-similar models which are asymptotic endstates for the general class [32 1 . 

Assuming a homothctic vector orthogonal to the cylinders of symmetry, we have obtained the standard form of 
the metric in cylindrically symmetric self-similar spacetimes. We have then applied this to the vacuum case and 
obtained solutions. In fact, the obtained solutions are all flat. We have explicitly shown that the spacetimes are part 
of the Minkowski spacetime with a regular or conically singular axis and with trivial or nontrivial topology. Although 
such spacetimes can emerge as the endstate of complete cylindrical gravitational collapse - for example, of cylindrical 
null dust (Tlj . this cannot be interpreted as evidence for the self-similarity hypothesis as the spacetimes are flat and 
therefore only trivially self-similar. 

Using the obtained self-similar expression for part of the Minkowski spacetime, we have argued that the C energy 
which is supposed to represent the gravitational energy per specific length of the cylindrically symmetric spacetime 
is subject to the choice of the translational Killing vector even if one chooses the same regular axis. We have also 
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(d) (e) (f) 

FIG. 3: The conformal diagrams of the two-parameter family of solutions for different values of k. The conformal diagrams of 
the two-parameter family of solutions for (a) < k 2 < 1/4, 1/4 < k 2 < 3/8, (b) n := 1/[2(1 - 2k 2 )} = 2, 3, 4, ■ • •, (c) n > 2 
but not an integer, (d) k < — l/v2, 1/V2 < k < 1, (e) k = 1, and (f) k > 1. The dashed lines denote cylindrical trapping 
horizons. The shaded and unshaded regions denote trapped and untrapped regions, respectively. The single circles denote 
timelike and spacelike infinities, while the double circles and double lines denote spacetime singularities. The solutions reduce 
to the Minkowski, the cylindrical Milne and the Kasner solutions for k — 0, 1/2 and —1/2, respectively. Choosing A ^ simply 
introduces a conically singular axis. For 3/8 < k 2 < 1/2 or n > 2, if n is an odd integer, there appears spacelike curvature 
singularity with r = [see (b)]; if n is an even integer, it is replaced by noncurvature but quasiregular singularity [see (b)]; if 
n is not an integer, then the spacetime admits not analytic but only a C' n ' extension beyond the null surface u = [see (c)], 
where [n] is the largest integer which is no greater than n. 
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discussed that there exists a cylindrical trapping horizon in the Minkowski spacetime and that the notion of trapping 
horizons might not be useful for defining black holes in cylindrically symmetric spacetimes - at least, not in the case 
where the marginal two-surfaces foliating the horizon are cylinders. 

Next, we have extended the analysis to the more general class of Einstein-Rosen waves, still respecting some kind 
of scaling behavior. Assuming a regular or conically singular axis, we have obtained a two-parameter family of non 
flat self-similar solutions, where the homothetic vector is not orthogonal to the cylinders in general. We have seen 
that the solution physically describes the interior of the exploding (imploding) shell of gravitational waves or the 
collapse (explosion) of gravitational waves depending on the parameter choice. Additionally, as a special case we 
have obtained a solution with a non azimuthal and non translational Killing vector which is not orthogonal to the 
cylinders. There is also a discrete subset of solutions which exhibit the collapse of gravitational waves developed from 
nonsingular initial data on a spacelike Cauchy surface. 

We conclude that self-similar Einstein-Rosen waves can describe nontrivial dynamics of gravitational waves if and 
only if the homothetic vector is not orthogonal to the cylinders of symmetry. Although the original proposal for the 
self-similarity hypothesis in general relativity is restricted in spherical symmetry, it is also likely that some of these 
self-similar solutions can describe the asymptotic behavior of more general solutions even in cylindrical symmetry. 
In fact, recent numerical simulations 27J strongly suggest that the asymptotic behavior of a dispersing gravitational 
wave within the null hypersurface t 2 = x 2 after the collapse of a dust cylinder is well approximated by a member of 
the family of solutions obtained here with k = —0.0206 and A ^ (see Figs. 8 and 9 of [27|). The present paper 
clarifies that this asymptotic solution belongs to the family of self-similar Einstein- Rosen waves with a non cylindrical 
homothetic vector and a conical singularity and that this asymptotic solution corresponds to gravitational waves 
inside the exploding shell of gravitational waves as shown in Fig. 3(a). Hence, it would be reasonable to generalize 
the self-similarity hypothesis - including in the context of gravitational collapse - as follows: under certain physical 
circumstances, solutions will naturally evolve to a self-similar form not only in spherical symmetry but also in a variety 
of symmetry classes. 
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APPENDIX A: TIMELIKE GEODESICS WITH INFINITE SPATIAL ACCELERATION IN THE PLANE 

SYMMETRIC CLOSED MILNE UNIVERSE 



The line element of the two-dimensional closed Milne universe is given by 

ds 2 = -dT 2 +T 2 d(p 2 , (Al) 

where < ip < 2tt and (p = and tp = 2n with the same T arc identified with each other. This is obtained 
from Eq. (|3.25|) or (|3 .31[) by omitting the two-dimensional plane part. The Milne universe is locally identical to the 
Minkowski spacetime. In fact, the coordinate transformation 

T = Tcoship and ( = Tsmh(p (A2) 

leads to the line element of standard form of the Minkowski spacetime 

ds 2 = -dr 2 + d( 2 (A3) 

In this (r, C)-coordinate system, the curves £ = and r = (/V are identified with each other, where V = tanh27r. 
More precisely, (t, 0) is identified with (r/yl — V 2 , tV/ y/l — V 2 ), where and hereafter we use the standard Cartesian 
coordinates (r, () of the Minkowski spacetime. See Fig. [4] 

We consider a past-directed timelike geodesic which starts from po ■ (t, £) = (tq , 0) with a unit tangent vector 
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FIG. 4: The closed Milne universe in the Minkowski spacetime. 



where < vo < 1 is assumed. This geodesic and the curve t = f/V intersect at the event 

VqTq Vv t 



Po ■ (r, C) 



V + v ' V + v a 



This event is identified with 



where 



pi : (r,C) = (ri,0), 



n := 



vg\/l- i'l 
V + v 



To 



(A5) 



(A6) 



(A7) 



We can determine the tangent to the geodesic at pi as follows. In coordinates (T,ip), the points po ■ (?b,0) and 
Pi : (T , 2tt) are identified. We must also identify the unit vector fields 

(A8) 



(A9) 



d 




d 




1 d 




i a 


df 


PO 


df 


pi 


flhp~ 


Po 


fd^ 



The coordinate transformation (|A2|) loads to 

d__rd_ <^d_ 1__3_ 

dT~Tdr + TdC T dip 
and the corresponding inverse relationship: 

d t d C d d 

Ifr = fdf ~ T*dip' d( = 



, d d 



(_d_ j_d_ 
fdf + T*dip' 



(A10) 



We can use this last equation to expand 



in terms of ^ and T invoke the identification 



revert 

to (t, C) coordinates by using (|A9|) and hence obtain the unit tangent to the geodesic at pi in coordinates (r, C) as 

1 d V\ d 



u 



(i) 



where 



vq + V 



vi 



1 + V vo ' 



(All) 



(A12) 
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It is also straightforward to show that in the transition from (tq,0) to (ti,0), a proper time of duration 



so ■= ri7~ To 

V + V 

elapses along the geodesic. 

Thus, this timelike geodesic goes through the points (r, £) = (ti, 0), (t2, 0), • • •. We can derive the recursion relations 



n+i 

Vi+l 
Si+l 



Ti, 



Here we introduce a new variable 



and rewrite Eqs. (|A13j) and (|A14|) as 



v + Vl 

v t + V 
1 + VV 

^y/r^yi-y 2 

Vi + 2V + v,V 2 1 



1 - Vi 



{l-Si)y/Si(2-5i) 

V + l-Si 
1- V 



1 + V(l - 5 t ) 



Prom Eq. (|A18[) . we have 



Si+i 



l-V 



S, i + V(l- St) 

Since < So < 1 and < V < 1, we have < <$i < <$o and hence < < 5;. This implies 

l-V 1* 

d . 



1 + V(l - S 



< s i+1 < 

Therefore, 

lim 5j = 0. 

i — >oo 

As for since < Si < So < 1 for i > 1, we have from Eq. (j A17|) 



Q < n+i = (1-8^(2 -Sj) < < L 



Therefore 



V + 1 - & 



lim Tj = 0. 

z — >oo 



V + l-5, 



(A13) 
(A14) 

(A15) 
(A16) 

(A17) 
(A18) 

(A19) 

(A20) 
(A21) 

(A22) 
(A23) 



Thus the three- velocity v of the timelikc geodesic which approaches p^ : (r, £) = (0, 0) becomes asymptotically the 
speed of light. The spatial part of the past-directed timelikc geodesic is infinitely accelerated in an approach to the 
origin : (r, C) = (0,0). 

We note however that p^ lies at a finite time in the past along the history of the geodesic. The total proper time 
that elapses along the geodesic is given by 



i=0 
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However, 

Si ~ v i + 2V + v i V 2 ' 

so the series converges: s < +00. 

Thus, there is no unique extension of the geodesies beyond (r, £) = (0, 0) and this behavior of geodesies is quite 
analogous to that around a conical singularity. This is solely due to the topological identification in ip and not related 
to the blow up of curvature. This corresponds to a quasiregular singularity defined by Ellis and Schmidt [30j . 



[1] A. Einstein and N. Rosen, J. Franklin Inst. 223, 43 (1937). 

[2] D. Kennefick, Travelling at the Speed of Thought. Einstein and the Quest for Gravitational Waves. (Princeton University 

Press, Princeton, 2007). 
[3] K. S. Thorne, Phys. Rev. 138, B251 (1965). 
[4] M. A. Melvin, Phys. Rev. 139, B225 (1965). 

[5] T. A. Apostolatos and K. S. Thorne, Phys. Rev. D46, 2435 (1992). 

[6] F. Echeverria, Phys. Rev. D47, 2271 (1993). 

[7] P. S. Letelier and A. Wang, Phys. Rev. D49, 5105 (1994). 

[8] T. Chiba, Prog. Theor. Phys. 95, 321 (1996). 

[9] A. Ashtekar, J. Bicak and B. G. Schmidt, Phys. Rev. D55, 687 (1997). 
[10] S. A. Hayward, Classical Quantum Gravity 17, 1749 (2000). 
[11] B. C. Nolan, Phys. Rev. D65 104006 (2002). 
[12] A. Wang, Phys. Rev. D68 064006 (2003). 

[13] K. Nakao and Y. Morisawa, Classical Quantum Gravity 21, 2101 (2004). 
[14] K. Nakao and Y. Morisawa, Phys. Rev. D71 124007 (2005). 
[15] Y. Kurita and K. Nakao, Phys. Rev. D73 064022 (2006). 

[16] K. Nakao, Y. Kurita, Y. Morisawa and T. Harada, Prog. Theor. Phys. 117, 75 (2007). 
[17] K. Nakao, D. Ida and Y. Kurita, Phys. Rev. D77 044021 (2008). 
[18] B. J. Carr, (unpublished). 

[19] B. J. Carr and A. A. Coley, Gen. Relativ. Grav. 37, 2165 (2005). 
[20] T. Harada and H. Maeda, Phys. Rev. D63 084022 (2001). 
[21] M. Snajdr, Classical Quantum Gravity 23, 3333 (2006). 

[22] C. Gundlach and J. M. Martm-Gar cia, Living Rev. Rel. 10, 5 (2007). URL (cited on 4th September 2008): 

http:/ /www.livingreviews.org/lrr-2007-5 
[23] J. Wainwright and G. F. R. Ellis Dynamical Systems in Cosmology, (Cambridge University Press, 1997) 
[24] M. Sharif and S. Aziz, Int. J. Mod. Phys. A20, 7579 (2005). 
[25] M. Sharif and S. Aziz, Int. J. Mod. Phys. D14, 1527 (2005). 

[26] L. V. Nolan, Cylindrically Symmetric Models of Gravitational Collapse, PhD thesis, Dublin City University, 2007 (unpub- 
lished). 

[27] K. Nakao, T. Harada, Y. Kurita and Y. Morisawa. rarXiv:0905.3"968l [Prog. Theor. Phys. (to be published)]. 

[28] H. Stephani, D. Kramer, M. A. H. MacCallum, C. Hoenselaers and E. Herlt, Exact solutions of Einstein's field equations, 

(Cambridge University Press, 2003). 
[29] J. M. M. Senovilla, Gen. Relativ. Grav. 30, 701 (1998). 
[30] G. F. R. Ellis and B. G. Schmidt, Gen. Relativ. Grav. 8, 915 (1977). 
[31] E. W. Hirschmann, A. Wang and Y. Wu, Classical Quantum Gravity 21, 1791 (2004). 
[32] C. G. Hewitt, J. Wainwright and M. Glaum, Classical Quantum Gravity 8, 1505 (1991). 



21 



ERRATUM 

We made a typographical error in the first term on the right hand side of Eq. (5.11). This is rectified as follows: 

ds = -[2(1 - 2/i 2 )p 2 e 2A y — - — J dmfo+f— — J ( v ^ _ J ^ + ( ~ 2~ J 

Also in the caption of Fig. 3, we erroneouly wrote that Fig. 3(b) showed the conformal diagram for n = 2, 3, 4, • • •. 
In reality, Fig. 3(b) shows the conformal diagram only for n = 3, 5, 7, ■ • ■. As for n = 2, 4, 6, • ■ -, the conformal diagram 
is given by the following figure, where unshaded regions denote untrapped regions. For this case, the extended region 
— v < u < is untrapped as is the original region < u < v. The null surface u = is a trapping horizon. The timclike 
surface r = is a regular or conically singular axis, while the spacelike surface r = is noncurvature quasiregular 
singularity. 
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